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Evidence for Fg Ground-State Symmetry of Cubic YbBi 2 
Probed by Linear Dichroism in Core-Level Photoemission 
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We have successfully observed linear dichroism in angle-resolved Yb^’'" 8^5/2 core-level photoe¬ 
mission spectra for YbBi2 in cubic symmetry. Its anisotropic 4 / charge distribution due to the 
crystal-field splitting is responsible for the linear dichroism, which has been verified by spectral 
simulations using ionic calculations with the full multiplet theory for a single-site Yb^^ ion in cu¬ 
bic symmetry. The observed linear dichroism as well as the polarization-dependent spectra in two 
different photoelectron directions for YbBi2 are quantitatively reproduced by theoretical analysis 
for the Ts ground state, indicating the Ts ground-state symmetry for the Yb®"*" ions mixed with the 
Yb^+ state. 


Rare-earth-based strongly correlated electron systems 
show various interesting phenomena such as competition 
between magnetism and unconventional superconductiv¬ 
ity, charge and/or multipole ordering, and the formation 
of a narrow (~meV) gap at low temperatures. AmoM 
them, YbBi 2 is known as a Kondo semiconductor [ll- 
[^, which has been recen^ recognized as a candidate 
for topological insulators p, as intensively discussed for 
another Kondo semiconductor, SmBe [S-Q- The mean 
valence of YbBi 2 has been estimated as ~2.9-|- by bulk- 
sensitive 3d core-level hard X-ray photoemission (HAX- 
PES) spectroscopy 0 - To discuss the mechanisms of the 
gap opening at low temperatures [1, flll - [l^ and the pos¬ 
sibility of a topological insulator, it is essential to verify 
the ground-state symmetry of the Yb®’*' [4/^® (one hole)] 
state determined by crystalline-electric-field (CEF) split¬ 
ting. Nevertheless, it is unclear for YbBi 2 . Its eightfold 
degenerated Yb®+ 4 / 7/2 levels are considered to split into 
two quartets [13 (normally, two doublets and a quartet) 
due to CEF in YbBi 2 . The ground state of the Yb®+ ions 
is commonly asserted to be in the so-called Fg symme¬ 
try but a possible accidentally de gene rated Fg-I-Fr 

ground state is not completely excluded [I 3 . In this Let¬ 
ter, we show evidence for the Fg symmetry of the Yb®’*' 
sites in the ground state of YbBi 2 on the basis of linear 


dichroism in angle-resolved core-level photoemission. 

Generally, it is difficult to experimentally determine 
the 4/ ground-state symmetry. Inelastic neutron scatter¬ 
ing is useful, but other excitations such as phonon exci¬ 
tations often hamper the observation of magnetic 4/—4/ 
excitations. Linear dichroism (LD) in 3d—to—4/ soft X- 
ray absorption spectroscopy (XAS) for single crystals is 
powerful owing to the dipole selection rules, as reported 
for Ce compounds [l9l - [^ . However, it is not applicable 
to compounds in cubic symmetry, in which there is no 
anisotropic axis relative to the electric field of the inci¬ 
dent light. On the other hand, the selection rules work 
also in the photoemission process while the excited elec¬ 
tron energy is much higher than that in the absorption. 
Furthermore, there is another controllable measurement 
parameter in photoemission, called the “photoelectron 
detection direction” relative to the single-crystalline axis, 
in addition to the polarization direction of the excitation 
light. Indeed, by using LD in 3d core-level HAXPES 
spectra, the Yb®+ 4/ ground state has been determined 
for tetragonal YbCu 2 Si 2 and YbRh 2 Si 2 [13 ■ TD in the 
core-level HAXPES for cubic Yb compounds is also ex¬ 
pected to be observed, as discussed below. 

In the case of Yb®+ ions in tetragonal symmetry, the 
eightfold degenerate J = 7/2 state splits into four dou- 
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FIG. 1 . (Color online) (a) Simulated polarization-dependent 
Sd^/2 photoemission spectra of Yb®"*" ions assuming the 
crystal-field-split ground state in cubic symmetry in the [100] 
direction, together with the corresponding experimental ge¬ 
ometry. (b) Same as (a) but for the photoelectron in the 
[ 111 ] direction. The 4 /-hole spatial distributions for the cor¬ 
responding states are also shown. 


blets as 

jri) =c|±5/2) + V1-c2|t3/2), (1) 

jr?) = -V1 -c2|±5/2) + c|t3/2), (2) 

|ri) = 6|±l/2) + vT^|T7/2), (3) 

\Tl) = ± 1/2)-b\T 7/2), (4) 

where the coefficients 0 < 6 <l, 0 <c<l defining the 
actual charge distributions and CEF splitting energies 
depend on the CEF parameters 5°, S 4 , i?!, Bg, and Bq 
in Stevens formalism [^ . In the cubic symmetry, the 
eightfold degenerate J = 7/2 state splits generally into 
two doublets and one quartet as 

|r6) = V^|± 7/2)+ 777121 T 1/2), (5) 

[Ey) = -73/2] ±5/2) +1/2] T 3/2), (6) 


I , r-777T2|±7/2) + 7^|Tl/2) .7. 

' 1 1/2] ± 5/2) + 73/2] T 3/2) ’ ^ ’ 

The Eg and E 7 states correspond to the tetragonal Eg 
(with b = 77 / 12 ) and E 7 (with c = 1/2) states, respec¬ 
tively. Since their 4/ charge distributions deviate from 
spherical symmetry owing to the CEF splitting even in 
the cubic symmetry, as shown in Fig. [ 1 ] it is natural 
to expect the observation of LD in core-level photoemis¬ 
sion for cubic Yb compounds. The Eg, E 7 , and Eg 4/ 
charge distributions are elongated along the [ 100 ], [ 111 ], 
and [110] directions, respectively. Actually, we have per¬ 
formed ionic calculations including the full inultiplet the¬ 
ory [ 2 ^ and the local CEF splitting using the XTLS 9.0 
program [ 2 ^. All atomic parameters such as the 4/-4/ 
and 3d-4/ Coulomb and exchange interactions (Slater in¬ 
tegrals) and the 3d and 4/ spin-orbit couplings have been 
obtained using Cowan’s code based on the Hartree- 
Fock method. The Slater integrals (spin-orbit couplings) 
are reduced to 88 % (98%) to fit the core-level photoemis¬ 
sion spectra M- 

Simulated polarization-dependent Yb^+ ^d ^/2 core¬ 
level HAXPES spectra in cubic symmetry at two photo¬ 
electron directions ([ 100 ] and [111]) are shown in Fig. [T] 
In the case of YbBi 2 , since the accidentally degenerated 
Eg + E 7 states could be a candidate for the ground state, 
the simulations assuming the Eg+E 7 ground state and its 
4/ hole spatial distribution are also shown. LD defined 
by the difference in spectral weight between the s- and p- 
polarization configurations is reversed between the [ 100 ] 
and [111] directions for all states displayed here. LD for 
the Eg ground state has the same tendency as that for the 
Eg ground state. On the other hand, LD for the Eg state 
is the smallest since the 4/ hole spatial distribution for 
the Eg state is the nearest to a spherical shape among 
these three eigenfunctions. LD assuming the Eg + E 7 
state is reversed to that for the Eg state. These simula¬ 
tions indicate that the symmetry of the Yb^’*' state can 
be determined by LD in the core-level photoemission. 

We have performed LD in HAXPES [ 2 ^, at 
BL19LXU of SPring -8 [1^ nsing a MBS Al-HE hemi¬ 
spherical photoelectron spectrometer. A Si(lll) donble- 
crystal monochromator selected linearly polarized 7.9 
keV radiation within the horizontal plane, which was far¬ 
ther monochromatized using a Si(620) channel-cut crys¬ 
tal. To switch the linear polarization of the excitation 
light from the horizontal direction to the vertical direc¬ 
tion, two single-crystalline ( 100 ) diamonds were used as 
a phase retarder placed downstream of the channel-cut 
crystal. The Pl (degree of linear polarization) of the 
polarization-switched X-ray after the phase retarder was 
estimated as —0.93, corresponding to the vertically lin¬ 
ear polarization component of 96.5%. Since the detec¬ 
tion direction of photoelectrons was set in the horizontal 
plane with an angle to incident photons of 60°, as shown 
in Figs, dja) and [Il[b), the experimental configuration 
at the horizontally (vertically) polarized light excitation 
corresponds to the p-polarization (s-polarization). The 
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FIG. 2. (Color online) Polarization-dependent Yb 3 ^ 5/2 core¬ 
level HAXPES raw spectra (solid lines) of YbBi 2 in the [100] 
direction and optimized Shirley-type backgrounds (see text), 
which we have subtracted from the raw spectra (dashed lines). 
The raw spectra have been normalized by the Yb^+ spec¬ 
tral weight. The same spectra in the expanded scale are also 
shown. 


excitation light was focused onto the samples using an el¬ 
lipsoidal Kirkpatrik-Baez mirror. To precisely detect LD 
in the Yb 3d core-level photoemission spectra, we opti¬ 
mized the photon flux so as to set comparable photoelec¬ 
tron count rates between the s- and p-polarization con¬ 
figurations. Single crystals of YbBi 2 synthesized by the 
traveling-solvent floating-zone method ii were frac¬ 
tured along the ( 100 ) plane in situ, where the base pres¬ 
sure was ~lxl0“^ Pa. The experimental geometry was 
controlled using a newly developed two-axis manipula¬ 
tor [^ . where the normal emission direction parallel to 
the [100] direction in Fig.[T][a) was changed to the photo¬ 
electron detection in the [111] direction in Fig.jljb) by az¬ 
imuthal rotation of 45° and polar rotation by ^55°. The 
sample and surface qualities were examined on the basis 
of the absence of any core-level spectral weight caused 
by possible impurities including oxygen and carbon. The 
energy resolution was set to 400 meV. The measuring 
temperature is 9 K, which is sufficiently lower than the 
excited state (> 250 K) [T^. 

The polarization-dependent Yb^"*" 8 ^ 5/2 HAXPES 
spectra in the [ 100 ] direction of YbBi 2 are shown in 
Fig. m A single peak at a binding energy of ~1519 eV 
and a multiple peak ranging from 1524 to 1534 eV exist 
in all the spectra. Since the 4/ subshell is fully occupied 
in the Yb^'*' sites with a spherically symmetric 4/ dis¬ 
tribution, the former single peak is ascribed to the Yb^+ 


states. The hnal states for the Yb^’*' components 

show an atomic-like multiplet-split peak structure in the 
1524— 1534 eV range. We show the same raw spectra in 
the expanded scale ranging from 1524.5 to 1530.5 eV in 
Fig.j^ A slight but intrinsic LD is seen in the highest and 
second-highest peaks in the raw spectra. The so-called 
Shirley-type backgrounds are also displayed in Fig. [51 We 
have optimized the backgrounds as follows: After the 
normalization of the background-subtracted spectra by 
both Yb^+ and Yb^’*' 3d^/2 spectral weights, the rela¬ 
tive Yb^+/Yb^+ contributions and the intensities in the 
high-binding-energy region of 1534 — 1540 eV become 
equivalent between the s- and p-polarization configura¬ 
tions. The reference binding energy on the higher side 
has been set to 1536.2 eV corresponding to the local min¬ 
imum of the raw spectral weight in the p-polarization. 
As a result, there are finite spectral weights at ~1536 eV 
in the background-subtracted spectra. However, these 
should be intrinsic owing to the overlap of the tails of the 
lifetime-broadened Yb^"*" 3 (i 5/2 main peaks and a plas- 
monic energy-loss structure at the higher binding energy. 
Note that we have confirmed the robustness of LD in the 
1524— 1531 eV region regardless of background intensity. 

A comparison of the polarization-dependent 
background-subtracted Yb^+ 8 ^ 5/2 HAXPES spec¬ 
tra of YbBi 2 and their LD with the photoelectron 
directions of [ 100 ] and [ 111 ] with the simulated ones for 
the Fg ground state is shown in Fig. [3| The highest peak 
is slightly stronger in the p-polarization configuration 
(p-pol.) than in the s-polarization one (s-poL), and 
the second highest peak is stronger in the s-pol. for 
both experimental and simulated spectra in the [ 100 ] 
direction. These tendencies are reversed in the data in 
the [ 111 ] direction. As shown in Figs.EKa) and[3][b), the 
observed LD and spectra are quantitatively reproduced 
by the simulations for the Fg ground state, for which 
the 4/ charge distribution is shown in Fig. [3][c). If the 
4/ ground state of YbBi 2 were in the Fg symmetry, 
LD would be much larger than the experimental one. 
The sign of LD for the Fy or accidentally degenerated 
Fg -f Fy ground state is completely inconsistent with our 
experimental results. Such a quantitative reproducibil¬ 
ity of the observed LD and core-level spectra by the 
simulations surely indicates the Yb^"*' (4/^^) ions in the 
Fg symmetry mixed with a small quantity (^ 10%) of 
the Yb^'*' (4/^^) component [l^ in the ground state of 
YbBi2. 

Our finding of the Fg 4/ symmetry for the Yb^’*' sites 
in the ground state is consistent with the prediction by 
a band-structure (local density approximation -|- on-site 
Coulomb repulsion, LDA-|-?7) calculatioin where the Fg 
and Fy states are on the occupied side Q. One might 
consider that the Fg and r 7 states are possibly mixed in 
the ground state due to the hybridization between the 4/ 
and valence-band orbitals at low temperatures well be¬ 
low the Kondo temperature (^ 240 K for YbBi 2 [IHJ)- 
As widely recognized for the local electronic structures 
of transition metal oxides, on the other hand, note that 
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the so-called CEF splitting actually seen in realistic ma¬ 
terials is a consequence of the anisotropic hybridization 
effects in addition to static ligand potentials. The latter 
would be much smaller than the former for YbBi 2 since 
the ligand-field potential on the Yb sites is to some ex¬ 
tent nearer to a spherically s ym metric one, being caused 
by the crystal structure [I|, y, 1^ where the Yb ion is sur¬ 
rounded by the truncated octahedron made of 24 boron 
ions as shown in Fig. [Sjd). The 4/-hole spatial distri¬ 
butions are elongated along the centers of the truncated 
octahedron faces for the Te and Tr states, whereas they 
are elongated along the edges of the hexagonal faces for 
the Tg state, which leads to the conclusion that the 4/ 
holes in the Fg state are relatively stabilized by the hy¬ 
bridizations compared with those in the Tg and F 7 states 
with different symmetries, as suggested by the LDA-|-t7 
calculation. Then, the possible Fg and/or F 7 state mix¬ 
ture would be experimentally insignificant in our data. 

In summary, we have successfully determined the 4/ 
symmetry of the Yb^+ sites in the ground state for cubic 
YbBi 2 as the Fg symmetry by LD in the Yb^"*" core¬ 
level HAXPES in two different photoelectron directions. 
Our result also suggests that the Yb^+ ion model un¬ 
der the effective CEF, in which the hybridization effects 
are implicitly taken into account, is suitable even for the 
valence-fluctuating system at low temperatures well be¬ 
low the Kondo temperature. The applicability of LD in 
the core-level HAXPES even to the system in the cubic 
symmetry (not restricted to systems with lower symme¬ 
try) demonstrated here would be promising for revealing 
the strongly correlated orbital symmetry of the ground 
state in a partially filled subshell, the charge distribution 
of which deviates from the spherical symmetry. 
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